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Abstract 

We do not present any original or new material. This is a tutorial addressed 
to students who need to study the microscopic derivation of the quantum- 
mechanical master equation encountered in many practical physical situations. 
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1 Introduction 

1.1 Motivation 

The term "master equation" (ME) used in this work means an equation of motion 

for the reduced density operator pA(t) of subsystem A which interacts with (usually 

much larger) another subsystem B. Experiences gathered during quite a few years 

of lecturing indicate that students find it difficult to understand the derivation and 

consequences of ME. The literature sources (at least those known to us) are usually 

quite brief and not easy to follow. 

The aim of this paper is to give a full presentation of the so-called microscopic 

derivation of ME together with detailed discussion of the underlying assumptions 

and approximations. Hence, this work does not bring any original material. It is 

just a tutorial which, hopefully, will help the students to obtain a better grasp of 

the ideas and concepts leading to one of important theoretical methods used in a 

variety of physical problems. 

The literature of the subject is quite rich. Entering the term "master equation" 

into Google Scholar returns about 30 000 links. It seems virtually impossible to give 

an extensive bibliography. Therefore, we concentrate only on several works which 

were essential in the preparation of this paper. 

1.2 Outline of the problem 

In order to present the ME technique, we first give the basic ideas which help, so 
to speak, to set the scene for further developments. We assume that the reader is 
familiar with fundamental concepts of quantum mechanics (given, for example, in 
the main chapters of the first volume of the excellent book by Cohen- Tannoudji et 



To start with, we recall that the state of the quantum system is, in the majority 
of practical applications, given by the density operator. This concept is introduced 
and discussed in virtually all handbooks on quantum mechanics, hence we mention 
only the basic facts which will be needed here. Namely, density operator for any 
physical system must have three essential properties 



The last inequality means that ( if) \ p \ ip ) > for any state | tp ) from the Hilbert 
space of the states of the given physical system. These properties may be phrased 
in terms of the eigenvalues A& of the density operator. Namely, A& G R, A& = 1, 
and Afc > 0. 

For a closed system with hamiltonian H, the evolution of the corresponding 
density operator is given by the von Neumann equation 



al&). 



Tv{p} = 1, 
P>0, 



P = P, 



hermiticity; 
normalization; 
semi — positivity. 



(la) 
(lb) 
(lc) 



ih- p(t) = [H, p(t)] 



(2) 
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which has a well-known solution 

p(t) = U(t,t )p(t ) rf(t,to), with U(t,t ) = exp(-^H(t-t )j , (3) 

where p(to) is a suitably chosen initial condition. Such an evolution is unitary pQ 

and obviously preserves all the necessary properties ([T]) of the density operator. In 

such a case everything is conceptually clear although necessary calculation may be 

quite involved or even requiring some approximate computational methods. 

The problem arises when we deal with a bipartite system A+B, consisting of two 

interacting subsystems A and B. Let us briefly outline the physical situation. We 

assume that the whole system A + B is closed and the total hamiltonian is written 

as 

H AB = H + V A b, where H = H A ®1 B + 1 A ® H B , (4) 

where Ha and H B are free, independent hamiltonians of each of the subsystems 
A and B. Vab is the hamiltonian describing the interaction between two parts. 
Let p AB (t) denote the density operator of the total system. Then, pAB^t) evolves 
according to the von Neumann equation (identical as ([2])) 

ift ^ PAB{t) = [ H A B, PAB{t) ] , (5) 

Now, one may ask, what is the problem? The point is that we are, in fact, inter- 
ested only in the subsystem A. Subsystem £>, for this reason or other, is considered 
irrelevant, although the A-B interaction certainly affects the evolution of A. More- 
over, in many practical cases, subsystem B is much larger (with many more degrees 
of freedom) and virtually inaccessible to direct measurements. Frequently, B is a 
reservoir [2, 3j and plays the role of environment upon which we have neither control 
nor influence. This may be very important in the context of quantum information 
theory [H [5] when the relevant subsystem is disturbed by the surroundings. More- 
over, problems of decoherence and irreversibility are intrinsically connected with the 
effects occurring in a subsystem influenced by an external reservoir (see [HE]). We 

shall not discuss these problems but focus on master equation technique. 

In the view of these brief remarks the question is, how to extract useful infor- 
mation on A from the general von Neumann equation ([5]). We stress that we are 
interested only in system A, so we need to find the reduced density operator 

PA (t) = Tr B {p AB (t)}. (6) 

The aim is, therefore, twofold: 

(i) extract the evolution of p A (t) from Eq. (|5]) for the entire system A + B. 

(ii) do it in a way which guarantees that properties (pQ) of p A (t) (as of any density 
operator) are preserved for any moment of time. 
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The solution to the stated problems is found in the so-called master equation 
technique. There are two, conceptually different but complementary, approaches to 
ME. 

The first one uses mathematically rigorous methods. Such an approach is pre- 
sented, for example, in Refs. [21 [3] (see also the references given in these books). 
Rigorous mathematics is obviously very important from the fundamental point of 
view. Mathematical theorems prove that ME for the reduced density operator of 
subsystem A follows from general von Neumann equation ([5]) and indeed preserves 
properties (UJ). In Refs. [21 0, SI it is shown that it is so, when ME attains the so- 
called standard (Lindblad-Gorini-Kossakowski-Sudarshan) form (Ref. [2], p. 8, Eq. 
(38); [3], Eq. (3.63) and 0], Eq. (78)) 

±p A (t) = ^[H, p A (t)] + Y^ajF iPA {t)F} - \[F}F U PA (t)] + \ (7) 

ij 

The first term leads the unitary (hamiltonian) evolution, while the second one is 

sometimes called a dissipator [3|. Operators constitute a basis in the space of 

operators for A subsystem (see [2112]). Finally, (X<i j IS cL positive definite hermitian 

matrix. Here, we do not go into the details of the derivation of the standard form 

([7]) of the master equation, we refer the reader to the mentioned references (see 

also [3]). Quite interesting derivation of standard ME is given by Preskill [5j. His 

presentation is somewhat heuristic but certainly worth reading, especially if one is 

interested in connection with quantum information theory, quantum channels etc. 

It is not, however, our aim to pursue formal mathematical issues. Our intentions are 

quite practical, so the reader may ask, why do we speak about the standard form of 

ME. The reason is as follows. 

The second approach to the stated problems is via "microscopic derivations". 

This occurs when we need to consider a specific physical situation when interacting 
systems A and B are known and well-defined. Then, we want to construct the cor- 
responding ME - equation of motion for the reduced density operator pA(t)- This 
is less formal and may be mathematically uncertain. If the microscopic derivation 
yields an equation in the standard form ([7]), we can say that the aim is achieved, 
because standard form ensures the preservation of the necessary properties of re- 
duced density operator pA(t)- Hence, both approaches are complementary. Formal 
but rigorous mathematical methods lead to standard form of ME which must be 

matched by equations obtained through microscopic derivation. 

All of the already mentioned references give (usually brief) account on the micro- 
scopic derivation of ME. These presentations seem to be difficult for students who are 
not acquainted with the subject and who seek the necessary introduction. Perhaps 
the most extensive microscopic derivation is given by Cohen- Tannoudji, Dupont-Roc 
and Grynberg [6j. This latter presentation is somewhat heuristic and, as it seems 
to us, leaves some nuances unexplained. There is, however, one more drawback. 
Namely, Cohen- Tannoudji et al do not compare their ME with the standard form. 
Therefore, essential question of positivity preservation remains untouched. 
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The derivation given here uses the concepts which can be found in the all ref- 
erences. Nevertheless, we are most indebted to Cohen- Tannoudji et al. Their ap- 
proach and especially their discussion of the time scales and employed approxima- 
tions strongly influenced our tutorial. We apologize if, at some places, we do not 
give proper references. Too many ones can distract the reader, which does not lessen 

our debt to all authors of cited literature. 

The scheme of this paper is summarized in the Contents, hence we feel no need 

to repeat it. We try to be as clear and as precise as possible. We focus our attention 

on the microscopic derivation of ME, leading to the standard form. Some issues are 

postponed to Auxiliary section, so that the main flow of derivation is not broken by 

additional comments, which can safely be given later. We hope, that the students 

who need a primer in the subject, will find our work useful and informative. The 

readers are invited to comment, so that the next version (if such a need arises) will 

be improved and, more readable. 



2 Evolution of the reduced density operator 
2.1 Introductory remarks 

We consider a physical system which consists of two parts A and B. We are interested 

only in what happens in part A which is usually much smaller than part B. The 

latter one we will call a reservoir (environment). We assume that the entire system 

A + B is closed. Then, its hamiltonian is specified as in Eq. ^ Some additional 

assumptions concerning both subsystems will be introduced when necessary. 

A previously, let pab {t) denote the density operator of the whole system A + B. 

The evolution of this operator is governed by von Neumann equation ([5]). Our main 

aim is to find the corresponding equation of motion for the reduced density operator 

PA = Tt{pab} for the subsystem A. Our starting point is provided by von Neumann 

equation which, after the transformation to the interaction picture, reads 

J t QAB{t) = ^ [V AB (t), Q AB (t)], (8) 

where, we obviously denoted 

Q AB {t) = e iHot / h PAB (t) e -W\ V AB {t) = e iHot / h V A Be- iHot /\ (9) 

with Hq given in Eq. Reduction of the density operator (as in ([6])) is preserved 
in the interaction picture (see Auxiliary sections) 

Q A {t) = Tv B {QA B (t)}. (10) 

Formal integration of Eq. (JSj) yields the following expression 

t+At 

g AB (t + At) = g AB (t) + i J dti [V AB (h),g A B(ti)], (11) 
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which gives the density operator at a later moment t + At, while the initial one at 
a moment t is assumed to be known. Iterating further and denoting 

AftAB(t) = QAB(t + At) - QAB(t), (12) 
we obtain, (similarly as in Ref.|6j) 

t+At 

Ag AB (t) = (J-) J dti [Vks(il), QABit)] 
t 

2 t+At 

+ (J-) J dt x J 1 dt 2 [V AB (tl), [V A B(t 2 ), QAB{t)} 



\ 3 t+At 

+ (jk) / dh ^ i* ^ [^AB(tl), [V A B(t 2 ), [V AB {tz) QAB{h)]. (13) 
t 

Higher order iterations will contain fourfold, etc., integrals and commutators. Let 
us note that in the last term we have time ordering t + At > t% > t 2 > > t. The 
above equation is rigorous, no approximations have been made. 

2.2 Weak-coupling approximation 

Weak-coupling approximation (discussed in all Refs.[2]-[6]), consists in retaining the 
terms up to the second order in interaction hamiltonian. Higher order terms are then 
neglected. Thus, we remain with 

t+At 

Ag AB (t) = (J-) J dtx [VAufa), Q A B{t)] 
t 

2 t+At t 

+ {jh) l dh J t 1(lt2 [VABfo), QAB(t)]- (14) 

t 

Alternatively, we can say that the obtained equation is valid in the second-order 

perturbation theory. Such an approximation requires a justification. This will be 

presented in the Auxiliary sections, now we focus on further steps of the derivation. 
Reduction of the operator qab (t) in the left hand side of (fTlj) poses no difficulties. 

Tracing over the reservoir variables (subsystem B) we obtain 
t+At 

Ag A (t) = (J^J J dh Tr B [V A B(ti), QA B (t)] 
t 

2 t+At t 

+ (J-) J dh£ 1 dt 2 Tr B [V A B(h), [V AB (t 2 ), Q AB {t)]. (15) 
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This expression has certain drawback. The the commutators in the right hand side 

contain full density operator @AB(t), and not the interesting (relevant) reduced one 

QA(t). To proceed, we need some more assumptions and approximations. 

One more remark seems to be in place. Subsequent iterations leading to Eq. (fT3l) 

are rigorous. In equation (]15p - which is approximate - there occurs the operator 

PAB(t), taken at the initial moment. The last term in the exact equation (|13p 

contains qab f° r moments earlier than the current moment t + At, but later than 

the initial instant t. This means that we neglect the influence of the "history" on 

the present moment. We shall return to the discussion of this point. 

2.3 Neglecting the initial correlations 

The key role in our consideration is played by the assumption that there are two 
distinct time scales [6]. The first one is specified by tb - typical time during which 
the internal correlations in the reservoir B exist. This will discussed in more detail 
later. Here we will only say that tb is such a time, that when it elapses, the state 
of the reservoir is practically independent of its initial state. The second scale is 
provide by time Ta- It characterizes the evolution (changes) of the operator QA(t) 
which is (are) due to the interaction with reservoir, and which may be specified by 
the relation 

A~t ~ T A 6A{t) - (16) 

Time Ta may be called the characteristic relaxation time of subsystem A. Let us 
note that we are speaking about interaction - the interaction picture we employ is 
thus, particularly useful. We make no statements about the rate of the free evolution 
of pa (in the Schrodinger picture), which is governed by hamiltonian Ha- Usually, 
the characteristic times of free evolution (the times of the order of ta ~ (Ha) a/ ft) 
are typically much shorter that Ta describing the influence of the interaction between 
subsystems [6]. 

Now, we assume that the introduced time scales satisfy the requirement 

t b < At < T A . (17) 

We have a fast scale (small tb) determining the decay of correlations within the 
reservoir and the second - much slower - scale defined by relatively long relaxation 
time Ta, characterizing the interaction between the two parts of the entire physical 
system. This may be phrased differently. We have assumed that the interaction is 
weak. Let V denote the average "strength" of this interaction. Uncertainty principle 
states that 

VTa ~ h => Ta ~ pp. (18) 

The condition tb <C Ta implies that 
h Vtr 

r B «T A ~- ^«1- (19) 

V n 
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In other words, we can say that spectral widths are the reciprocals of characteristic 
times, so the condition tb <C Ta means that the spectral width of the reservoir 
energies must be much larger than the spectral width of the interaction between 
subsystem A with reservoir. Further discussion and justification of our approxima- 
tions is postponed to Auxiliary sections. Here we focus on the derivation of the 
master equation. 

The adopted assumption tb <C Ta allows us to make the approximation which is 
sometimes called the Born one [3]- [6]. Initial density operator for the whole system 
A + B can always be written as 

QAB{t) = QA(t) ® QB{t) + Qcorel(t), (20) 

where QA(t) and qb(P) are the reduced density operators for two subsystems. The 
state of the whole system consists of a factorizable part qa (t) <8> QB(t) and the en- 
tangled part Qcorei{t)-, describing the interaction-induced correlations between the 
subsystems. Equation (fT5j) gives us the change AgA(t) = QA(t + At) — ^4 (t ) , 
hence informs us about changes occurring in the time interval At. Assumption 
that tb <C At allows us to neglect the mentioned correlations. As previously, we 
postpone the discussion for Auxiliary sections. At present, we assume that 

QAB{t) « QA{t)®QB{t). (21) 

By assumption, the reservoir (environment) is large, its correlation time is very 
short, so the reservoir's relaxation is fast. We may say that before any significant 
changes occur in subsystem A, the reservoir would have enough time to reach ther- 
modynamic equilibrium. As it is known from statistical physics such state is given 
as 

°B = ^2p(z)\z)(z\ where p(z) = - exp (-jrj; J • ( 22 ) 

Z. 

The quantity Z is a partition sum 

Z=I>(-^). (23) 

States I z ) and energies E z are the eigenstates and eigenvalues of the reservoir hamil- 
tonian, which can be written as Hb = ^2 Z E z \z){z\. As a consequence we conclude 
that 

[a B , H B ] =0, (24) 

so we can say that operator gb is stationary - does not change in time. It is worth 
noting that we could have reversed the argument. First require stationarity, as ex- 
pressed by (|24l) which would entail relations (|22|) and (|23j) . Moreover, commutation 
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relation (|24p implies that the reduced density operator ob is of the same form both 
in Schrodinger and interaction pictures. 

Due to these remarks operator g B (t) appearing in Eq. ([2Tj) is simply replaced by 

ob- Therefore, in Eq. (fT5|) we make the replacement qab = QA(t) (g) ob- So, we now 
have 

t+At 

AfiU(*) = (J-) J dh Tr B [VAB(ti), QA{t)®a B ] 
t 

2 * +Ai 

+ (J-) J dhj 1 dt 2 Tr B [ V AB (ti), [ V AB {t 2 ), Q A {t) ®* B ]. (25) 



The employed simplification facilitates computation of the remaining traces. How- 
ever, to proceed effectively, we need to specify the interaction hamiltonian. 

3 Interaction hamiltonian and its properties 

3.1 The form of V AB (t) 

Our next assumption concerns the shape of the interaction hamiltonian. It is taken 
as (similarly as in the given references) 

Vab = Yl A *® X <* = A i® x w (26) 

a a 

where A a are operators which act in the space of the states of subsystem A, while 
operators X a correspond to space of the reservoir's states. Operators appearing in 
the definition (I26D need not be hermitian (each one separately). Only the hamilto- 
nian Vab must be hermitian. That is why we have written the second equality. We 
can say that to each nonhermitian term A a ® X a corresponds the term Aa (g) X a , 
and the latter appears in the sum Vab, but with another number. In Auxiliary 
sections we will argue that it is not any limitation. It is only important that the 

whole hamiltonian Vab must be hermitian. 

Operators A a and X a act in different spaces so they are independent and com- 
mute. In the interaction picture we immediately have 

V AB {t) = Mt)®x*(t) = 4(<)®^<L(t), (27) 

a a 

with 

A a (t) = e iHAt/h A a e- iHAt/h , X a (t) = e iHBt/h X a e' iHBt/h . (28) 

Rules of hermitian conjugation imply that the conjugate operators transform to 
interaction picture in the exactly the same manner as the initial ones. 
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We now make one more assumption about reservoir. We have already taken 
q b (t) a B . Here, we assume that in the Schrodinger picture 

(X a ) B = Ti B {X aPB (t)} = Tr B {X a a B } = 0. (29) 

This assumption easily transforms to interaction picture 

( X a (t) ) B = Tv B { e iHBt ' h X a e - iHBt ' h a B } 

= Tv B { X a e- iHBt ' h a B e lHBt ' h } = Tr B {X a a B } = 0, (30) 



which follows from commutation relation (|24p . cyclic property of trace and (|29p . 
This is rather a simplification, not a restrictive assumption. This will be clarified 
and explained in Auxiliary sections. Eq. ([29]) (leading to ([30]) ) allows us to see that 
the first term in the ME (|25p is, in fact, zero. Indeed 



Tr B [VAB(h), 0A(t)®a B ) = TV B A a (t x ) ® X a (t x ), Q A {t)®a B 

L OL 

= YX Mh)QA(t) Tr B [X a {h) a B ] 



QAA a (ti) Tr B 



a B X ( 



o. 



(31) 



Both traces are equal (cyclic property), nevertheless this expression need not be 
zero, because operators of the A system need not commute. If requirement (|29p 
is not fulfilled then the above average may not vanish. Assumption (|29p an d its 
consequence (|30p fortunately give zero, thus the first term of Eq. (|25p vanishes and 
we remain with the master equation 



1 

ih 



2 i-t+At ft! 

dt\ I dt2 Tt b 



V AB {ti), [V AB (t 2 ), Q A (t)®a B ] ■ (32) 



Expanding the commutators is simple. Moreover, one easily notices that there are 
two pairs of hermitian conjugates. Hence we have 



AgA(*) 
At 



1 



t+At 



h 2 At 



dU 



dt 2 Tr B { V AB {t2) (QA(t) ® a B ) V AB (h 



V AB {ti)V AB {t 2 ) [Q A (t)®0 B )} 



+ 



.C. 



(33) 



We can now use hamiltonian (|27p and perform further transformations. It can be, 
however, shown [2]-[l] that this equation does not guarantee that the positivity 
of the density operator QA(t) is preserved. It appears that the so-called secular 
approximation is necessary. To perform it effectively, it is convenient to present the 
interaction hamiltonian in a somewhat different form. 
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3.2 Operators A a (Cl) 

Let us write the hamiltonian of the subsystem A as 

H A = Y^^a\a){a\. (34) 

a 

States | a ) constitute the complete and orthonormal basis in the space of states of 
the subsystem A. The eigenfrequencies uj a may or may not be degenerate. We allow 
oj a = uj b for a 7^ b. At present it suffices that we distinguish different kets | a ) 
solely by their "quantum numbed a. Similarly as in Refs.pU 0], we now define the 
operators A a {Vt) via the following relation 

A a (n) = S(u ba - n) | a )( a \ A a \ b )( b |. (35) 
a,b 

This representation may be called the decomposition of operator A a into eigen- 
projectors of hamiltonian H A - Delta 5{oJba — ^) is of the Kronecker type, that is 

x , n ; / for u) ba / VL, . . 

^ " n) = { i for u ba = n, (36) 

In our considerations we allow for nonhermitian operators A a . Hence, definition 
(f35|) is augmented by the following one 

4 (ft) = S ( u ba -n)\b)(b\At t \a)(a\ 

a.b 

= Y,8(e>c*-to)\a)(a\At [ \b)(b\, (37) 

a.b 

because it is always allowed to interchange the summation indices a <-> b. We stress 
that A a (Q) contains Bohr frequency uj ba , while in Aa(Q) we have uj ab = —u>ba- The 
following relation seems to be quite obvious 

^<5(w fe ,-n) = l, (38) 
n 

since out of all fl's one will exactly match Ukn- As a consequence we obtain 

Y,A a (Q) = A a . (39) 
n 

Indeed, from definition ([55]) . relation ([55]) and due to completeness of states | a) we 
get 

£j4«(n) = J3J3a(w te -n)|a)(a|i4a|&)(&| 

f2 f2 a, 6 

= I a)(a I A» I 6){6| = l>l a l = i4 a . (40) 
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Relation f|39|) implies that the interaction hamiltonian (in Schrodinger picture) can 
be written as 

v AB = Y; A «® x « = Y.Y. A ^® Xa - ( 41 ) 

a Q a 

Similarly as above we can also show that 

= 4, (42) 

n 

and 

vab = = E E A «w ® x l ( 43 ) 

Q a 

Using definition (|35h we find the operator Aj(Cl) in the interaction picture 

A a (n) = e iHAt/h A a (Q) e- iHAt/h = e~ int A a {n), (44) 

because e lHAt ^ h \ a) = e %u)at \a). Linking expressions (I4ip and (I44p we write the 
interaction hamiltonian in the interaction picture 

v AB (t) = E E e ~ int ® 

n a 

= E E e ^ 4(^)®%) = (*) • ( 4 5) 

a 

Before starting to analyze ME (j33|) . let us notice that operators A a (Q) possess 
some interesting properties [3]. Discussion of these properties is moved to Auxiliary 
sections. 

3.3 Further analysis of master equation 

We return to master equation f|33j) . Interaction hamiltonian Vab^) is taken as in 
the first part of Eq. (f4"5j) . while Vab(£i) is represented by its hermitian conjugate 
according to second part of (|45p . This gives 

t+At 

Ag> A (t) 



At h 2 



<nt3 4a(fi) ® X^(t 2 ) ® a B ]e in ' h AUn') ® X+(t x ) 



- e m '*Mt(nO®^(ti)[e^^(n)®^(t 2 )]ex(t)®ffB} + e.c. (46) 

Computing tensor products we remember that partial trace is taken only with respect 
to reservoir variables. Moreover, we note that these traces are the same (cyclic 
property). Therefore we denote 

G a0 {t x -t 2 ) = Tr B { Xl (h) Xp{t 2 ) a B } . (47) 
S.K, J.C-K 11 
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Finally we rewrite the arguments of the exponentials as ifl 't\ — iQ,t 2 = i(P- ' — Q)t\ + 
itt(t! - t 2 ). Thus Eq.(0S]) becomes 

t+At 

{ Jt a ,p n,n> 
Ap(tt)Q A (t)Ai{£l') - Ai{Sl')A p {Q)Q A {t)\ + H.C. (48) 

The quantity G a p{t\ — t 2 ) is called the correlation function of the reservoir. We will 
briefly discuss its properties. 

3.4 Correlation functions C7 Q/ g 

Let us focus on the functions defined by the right hand side of Eq. (|47p . that is 

G a p{h,t 2 ) = Tr B {xl(h) X (t 2 ) a B } ■ (49) 

These are the functions of two variables and it is not, a priori, clear that they 
depend only of the difference r = t\ —t 2 . Before discussing this fact, let us note that 
G* g(ti, t 2 ) = G/3 a (t 2 ,ti). To prove this relation, we recall that Tr^{^4} = Tr B {A' }, 
so that the definition (149p gives 



G* af3 (h,t 2 ) = Tr B [a B xl(t 2 ) X a (h) } 

= Tr B {x\{t 2 ) X a {h) a B } = Gp a (t 2 ,h). (50) 

Now, we will show that the function G Q p(ti,t 2 ) is indeed a function of the dif- 
ference of its arguments. The key role plays the fact that the state of the reservoir 
(density operator a B ) is stationary (does not change in time). Explicitly using the 
interaction picture we get 

G a/3 (ti,t 2 ) = Tt B { {e iHBtlth Xl e - iH ^/ K ) (e iHBt ^ h Xp e ~ iHBt2 / H ) a B } . (51) 

The trace is cyclic and a B commutes with hamiltonianem H B so we conclude that 

G a/3 (h,t 2 ) = Tr B { j"*^^'* X^ e -iH B {H-t 2 )/n Xp - g j 

= Tr B {xl(t 1 -t 2 )X /3 (0)a B } = G a(3 (r = h - t 2 ), (52) 

for two moments of time t\ > t 2 . Reservoir's correlation function effectively depends 
only on one variable. This fact is denoted by a bar over the symbol of correlation 
function. Thus we write 

G aP {t x ,t 2 ) = G a p{r) = Tr B [xUT)X p a B } (53) 
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Such correlation functions are called stationary. Sometimes the concept of station- 
arity means invariance with respect to time translation. Indeed, for arbitrary time 
T we have 



G a p{h+T,t 2 +T) = ^((ti + Tj-Ctz + T)) = G a/3 (ti-i 2 ). 



(54) 



This property of the correlation functions is a straightforward consequence of the 
stationarity of reservoir's density operator erg. 



4 Discussion of times 

Preparing this section of our tutorial we greatly benefited from the analogous dis- 
cussion by Cohen- Tannoudji et al [6j. To large extent we follow their reasoning, 
trying to elucidate some less obvious points. Admitting this, we will refrain from 
giving multiple references to their work. 

4.1 Limits of the integrals and Markov approximation 

In master equation (I48p one integrates over the triangle ABC which is shown in 



Fig. [TJ Firstly, one computes the integral over dt 2 in the range from t to t±. This 
is indicated by thin vertical lines (at left). Next, one sums such contributions by 
integrating over dt\ from t to t + At. The integrand in (|48|) contains correlation 
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Figure 1: Left figure presents the integration region in the double 
integral in Eq. ([4"5|) . Right figure illustrates the change to new 
variables r = t\ — ti and t\. Other explanations are to be found 
in the text. 



functions of the reservoir which depend on the difference r = t\ — ti- We stress that 

we always have t\ > so that r > 0. The integration over the triangle can be 

performed in another manner. 

Let us consider the geometry (right graph in Fig. [1]). Along the diagonal AC we 

have t\ = t 2 , so r = t\ — t 2 = 0. The straight line l\ has (in t\ and t 2 variables) the 
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equation t 2 = t\ — r, where r is fixed, since (— r) is the coordinate £2 of the point 
where the discussed line intersects the axis t 2 - Then, for the line 1^ (passing through 
point B) r is also fixed (by the same argument, as in the case of line On Z^, at 
the point B we have t\ = t + At and t 2 = t. Thus, we have r = At. Parametr r 
specifies the skew straight lines (parallel to the diagonal AC) and passing through 
triangle ABC. Integration over the triangle ABC is now done as follows. We fix 



) 

t + At 


I; 


t 2 = h 
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Figure 2: Illustration of the change of integration variables in 
equation (|48|) . Transformation to new variables r — t\ — ti \t\. 
Other explanations - in the text. 



r G (0, At) and we move along the segment A'C (see Figj2]). Variable ti runs in the 
interval from t + r to t + At. So, we first integrate over dt\ from t + r to t + At (along 
the segment A'C). Next, we integrate over dr from zero to At. In this manner we 
sum the contributions from all skew segments covering the triangle ABC. Therefore, 
we can write 

t+At t+At 
rti rAt r 

dti dt 2 = dr dti, (55) 

Jt Jo J 

t t+T 

while we remember that r = tj — t 2 (or t 2 = t\ — r). Performing the discussed 
changes of integration variables in Eq. (|48p . we get 

At t+At 



At h 2 At 



[dr [ d tl J2J2 * i{Q '~ m ^ T G a ,(r) 
t+r n.«' 

Ap(n)e A (t)Ai(n') - AUn')A^n)g A (t)} + e.c, (56) 



We recall that the considered time intervals satisfy the requirement At ^> tb (which 
will be discussed in detail later). If it is true, then the main contribution to the 
integral over dr in Eq. (|56p will come from the region in the neighborhood of < 
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r < tb <C At. Geometrically, this corresponds to a narrow belt which is parallel 
to the diagonal AC and lies just below it. It follows, that outside this region the 
reservoir's correlation functions practically vanish (decay to zero). Therefore, we 
will not make any serious error moving the upper limit of integration over dr to 
infinity. Moreover, since only small r's contribute significantly, the lower limit of 
the integral over dt\ may be approximated simply by t, so only a small "initial" 
error might be introduced. With these approximations equation (|56p yields 

t+At 

Ag A {t) 



n,C a,/3 ^° { 

A p (Q)g A (t)Al(Q') - AUn')Ap(n)Q A (t)] + H.C. (57) 



At /i 2 _ 

n,n' a,/3 



Introducing the quantities 

t+At 



J(n'-n) = — J dh exp[i(Q' -Sl)ti], (58a) 

W a/3 (Q) = dr e tnr G aP {T) = dr e^ T Ti B {xi{T)X p a B }, (58b) 
Jo Jo 

we rewrite Eq. (|57p as follows 

x [A p (n)e A (t)Ai(n') - Al(n')A p (n)Q A (tj\ + e.c. (59) 



This equation specifies the rate of change of the reduced density operator g A (t) 
within the time interval (t,t + At). The quotient Ag A (t)/At can be treated as an 
averaging 

Ag A (t) _ g A {t + At)-g A {t) _ 1 f t+At ^ dg A (h) (6Q) 



At At At J t dh 

This averaging results in smoothing all very rapid changes of g A (t) which may occur 
during the interval (t, t + At). In principle we should account for such rapid changes. 
We do not do that because right hand side of Eq. (l59p contains g A (t), while the left 
hand side represents the smoothed rate of change. This rate depends on the density 
operator g A in the past, that is at the moment when the smoothed evolution was 
started. So our next approximation consists in replacing the smoothed rate by a 
usual derivative. In other words, the variation at an instant t (that is the derivative 
dg A {t) / dt) is connected with the value of g A (t) at the very same instant. This 
approximation allows us to use a usual derivative at the left hand side of (|59p . This 
approximation sometimes is called a Markovian one since it connects the variations 
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of some physical quantity with its value at the same instant, independently from 

the values which this quantity had at earlier moments. We can say that Markovian 

approximation consists in neglecting the influence of the history of the physical 

system on its current state which fully determines the presently occurring changes. 

In some literature sources this approximation is also called the coarse-graining one, 

because small and rapid fluctuations are neglected when the evolution is investigated 

on a much longer time scale specified by At. 

With all the discussed approximation our master equation (|59[) becomes 

" Q,n' a,/3 

x [A p (Q)g A (t)Al(n') - 4(fi')4g(n)ftA(t)] + H.C. (61) 
4.2 Schrodinger picture 

At this stage, we return to the Schrodinger picture and we insert 

g A {t) = e iHAt / h PA {t) e- iHAt l h , (62) 



into Eq. (|6ip . When computing the derivative at the left hand side we reproduce 
the free evolution term. Thus, from (|61j) we get 



e iH A t/h(d\ e _ iHAt/h = _ ^ iH A t/n^ H ^ pA{i) ^ e - iHA t/n 

\ Liu J I v 

+ UEE J ( Q '- n)W^(fi) [ A P (Q) e iHAt ' h p A (t) e- iHAt ' h Al (Q') 
^ n,n' a,/3 

- Al(n')Ap(n)e iHAt / h p A (t) e~ iHAt ^ + M.C |. (63) 

Multiplying on the left by e ~ lHAt / h and on the right by e lHAt / h , we use relation (|4"4"|) 
and its hermitian conjugate (for negative times). This yields 



d PA (t) =-i[H A , PA (t)]+ I^Y,Y, J ^'-^ w ^ e 



dt h 



i(Q-Q')* 

A p (n)p A (t)Al(n') - Ai(n')Ap(n) PA (t)] + e.cj. (64) 

4.3 Integral J(0' — O) and secular approximation 

Our master equation contains the integral J(f2' — O) defined in (|58aj) . Its computa- 
tion is straightforward. Denoting temporarily x = f2' — fi, we get 



t+At ixti 



./(.,•) / dt, = e -t+-A t /2 !^V^Z = e *rt F(x)) (65) 



(^) 
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where we have introduced (as in [6]) a function specified as 



F(x) 



JxAt/2 



sin 



2 ) 



Due to the obtained results we can write 

j(o'-o) = e i{n '- n)t F(n' -n). 



(66) 



(67) 



Inserting the computed integral into (|64p we note that the exponential factor cancels 
out. Hence 



PA(t) 



[h a , PA (t)]^±Y,Y, F ( Q '- n ) w °pW 



A /9 (n)p A (t)4(n) - Ai(p) a p (q) PA {t) 



+ 



The sense of function F(£l' — Q) must be carefully considered. We proceed along 
the lines similar to those in [6]. It is easy to see that function |-F(x)| has a sharp 
maximum for x = 0' — f2 = 0, where it is equal to unity. 




Figure 3: The graph of the modulus of the function F(il ' — f2) which appears in 
(|68[) . If the time At is sufficiently large then the graph has a sharp maximum for 

ct' = a 



Zeroes of this function correspond to 

1 a 2n7r 

T^xAt = mr x = — — . (69) 

2 At y J 

If the time At is sufficiently long then the central maximum is very narrow. The 
question is what does it mean "sufficiently long time". Let us consider two possibil- 
ities. 

1. If x = \Q' — £l\ <C (At) -1 , the argument of function is very close to 

zero, its value being practically one. 
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2. If x = — Q,\ 3> (At) 1 (Bohr frequencies are significantly different) then 
|F(x)| is close to zero, as it is seen in Fig El 

We conclude that the terms at the right hand side of master equation (I68p containing 

the operator products Aa(Q, ') Ap(Q), for which |f2' — Q,\ 3> (At) -1 practically do 

not contribute to the evolution of the density operator pA(t). According to the 

first possibility above, significant contributions come only from such couplings that 

operators Aa(£l') and Ap{Q) have practically equal Bohr frequencies. 

As we know, time Ta is a characteristic relaxation time in subsystem A due to 

interaction with reservoir and it satisfies the estimate At <C Ta (we discuss it later). 

It can be argued (see also [6]) that the terms in master equation (|68p . in which 

|0/ — O) ~ (At) -1 also give very small contributions, so that they can be neglected. 

As a result of all these approximations, we may say that only those terms in right 

hand side of master equation (I68j) contribute significantly for which fi' — 0,1 = 0. 

Such an approximation is called the secular one. It allows us to replace the function 

F(£l' — O) by the Kronecker delta 5(0,' — f2) defined as in (I36p . It reminds us that 

only the terms satisfying the requirement (O' = S7) give nonzero contribution. All 

these arguments lead to master equation of the form 



d , . i 

5 = - S 



[ Ha, PA(t) ] + (bEE*'" ) w «eW 

> ™ fl O ' o 



a,/3 fi'.n 

A p (Q) PA (t) At(O') - 4(fi')4,(«)^(i)] + H.c|. (70) 

The presence of the discussed Kronecker delta simplifies one of the summations, 
which gives 

j PA {t) =- l -[H A ,PA{t)]+ I ±J2Y, w <*pW[Mn) P A(t)AUn) 

- Ai(n)A (n)p A (t)] +M.C |. (71) 

The fundamental part of the microscopic derivation of the master equation is fin- 
ished. We shall perform some transformations which have important, but rather 
cosmetic character. We want to transform master equation into the so-called stan- 
dard form. All other comments mentioned many times, left to Auxiliary 
sections. 



5 Standard form 
5.1 Introduction 

Standard form is important, because it can be shown (see [H O 0]) that this form 
guarantees preservation of hermiticity, normalization and, first of all, the positivity 
of the reduced density operator pa- If our master equation (fTTj) can be brought 
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into the standard form, then we can be sure that all the necessary properties of the 
density operator of subsystem A are indeed preserved. Obviously, the first term in 
the right hand side of equation (|7ip describes the unitary evolution, hence we shall 
concentrate only on the second term. Writing explicitly the hermitian conjugates, 
we have 



Q a,/3 



+ 



^EE^^)f^(°)^(*)4^) " PA(t)Al(n)A a (n)], (72) 



n a,/3 



because operator pA(t) is hermitian (the proof that hermiticity is preserved will be 
presented in Auxiliary sections). In the second term we interchange the summation 
indices a <-> (3 which gives 



d_ 

dt 



PA(t) 



1 



n a,(3 



+ M^)PA(t)Atm - p A (t)AUn)Ap(n) 



h 2 

fl a,/3 

5.2 New notation 

For further convenience we introduce the following notation 



(73) 



r a0 (n) = w aP (n) + w£ a (n), 



(74a) 
(74b) 



The matrix r Q/ g(Q) is hermitian and positively defined. The latter property is dif- 
ficult to prove. It requires some advanced mathematics and we take this fact for 
granted. The readers are referred to literature [21 [3|. Hermiticity of T a p follows 
directly from the definition (|74ap . Indeed, we have 



r^(n) = w*,(n) + w Pa (n) = r^o). 

The second matrix A a ^(Q) is also hermitian. From (|74b|) it follows that 



(75) 



A* a/3 (Q) = - 1 [W^(fi) - W Pa (Q)] = | [W Pa (Q) - W*p(n)] = Ap a (n). (76) 

Let us focus on the method of computation of elements r Q/ g(0). As it will be 
shown, elements A Q/ g(r2) are less important. To find T a p we need quantities W% a . 
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Conjugating definition (|58b|) we find that 



dr e lilT Ti B {Xl(T)X a a B } 

POO ____ 

/ dr e~ inr Tr B {xlXp(r)a B } 
Jo 

POO 

/ dr e~ inr Ti B {e- iHBT / n Xle iHBT / n X^a B }, 
Jo 



(77) 



where we used relations ([24"|h ([28]) and cyclic property of trace. Changing the inte- 
gration variable t — > — r, we have 



dre^ Tr B {xl(T)Xpa B }, 



(78) 



The integrand is identical as in the definition (|58bp . only the integration limits are 
different. Combining this with (I58bl) . we get 



/OO POO 
dr e inr Tr B {Xi{r)Xpa B } = / dr e inr G a(3 (r). 
-OO J — OO 



(79) 



The elements r Q( g(r2) are the Fourier transforms of the corresponding correlation 

function of the reservoir. 

Matrix A a p(£l) does not have such a simple representation. From the definition 

(|74bj) and the second relation in (|77l) 



(80) 



A Q/3 (0) = i [ / dr e inT Tr B {xl (r)X p a B } 



POO ____ 

/ dr e~ inr Tr B {XlX p (r)a B } 
Jo 



5.3 Standard form 



Inverting relations (|74p we express elements W a p via and A a g. After simple 
regrouping of the terms in Eq. (|73p we find 



d_ 

dt 



PA(t) 



= ^EE r aff {n){A p (Q) PA (t)A^(Q) 

4(0)^(0), paW 



-iA a/3 (0) [4(0)^(0), ^(t)] }. (81) 
Let us note that the last term is a commutator, so we define 

= \ EE 4^)4(0)^(0). (82) 



a,P 



Taking into account hermiticity of matrix A a g(u) (changing the names of the sum- 
mation indices when necessary) we can easily show that the operator H^s is a ls° 
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hermitian. Returning to full master equation, that is to Eq. (|71|) . we conclude that 
the term containing Hls in (|8ip can be connected with the free hamiltonian one. 
In this manner, we finally have 

j t PA(t) = -^[H A + H LS , p A (t)] 

+ ^EE r^(n){^(j2) P x(t)4(n) 

n a,/3 

-i[4(n)^(n), PA (t)) + }, (83) 

which coincides exactly with the standard form of the evolution equation for the 
reduced density operator pA(t) which describes the state of the subsystem A inter- 
acting with reservoir B. This allows us to be sure that hermiticity, normalization 
and positivity of the operator pA(t) are indeed ensured. Finally, let us remark that 
operator Hls which gives a contribution to the hamiltonian (unitary) evolution, usu- 
ally produces small shifts of the eigenenergies of the subsystem A. That is why, in 
many practical applications, this term is simply omitted. This explains our previous 
remark that matrix A Q/ g is less important than T a p. Obviously one can construct 
operator Hls and investigate its influence on the unperturbed energy levels of the 
subsystem A. Small energy shifts of eigenenergies of subsystem A are qualitatively 

similar to the well-known Lamb shifts, which clarifies the employed notation. 

The obtained master equation (|83p is an operator one. In practice, we frequently 

need an equation of motion for the matrix elements of the reduced density operator 

PA(t)- It seems to be natural to use the energy representation, that is to consider 

matrix elements of pa (t) calculated in the basis { | a ) } of the eigenstates of the free 

hamiltonian Ha (see Eq. (j34p ). This will be done in the next section. 

5.4 Energy representation 

When analyzing master equation in the basis of the eigenstates of free hamiltonian 
we need to be careful. The reason is that the commutator in (|83p contains an 
additional term, namely the Lamb-shift hamiltonian. One may argue that this 
changes the hamiltonian and a new basis should be found (a basis in which Ha+Hls 
is diagonal). We will, however, proceed in the spirit of the perturbative approach. 
We will treat Hls as a small perturbation which, at most, will yield small energy 
shifts. Therefore, the set {| a)} of eigenstates of the unperturbed hamiltonian Ha 
can be used as complete and orthonormal basis. Working within this scheme, we 
can easily construct master equation (equation of motion) for matrix elements of 
the density operator for subsystem A. We will suppress the index A since it should 
lead to no misunderstanding. Taking matrix elements p a b(t) = ( a \ PA(t) \b) and 
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expanding the anticommutator term we obtain 
j t Pab(t)=- l -{a\[H A + H LS , p(t)]\b) 

+ ^EE r ^){H^)p(*)^)i&> 

- \ <a|4(Q) A P (Q) P (t)\ b) - \ (a\p(t) 4(0) Ap(fl)\ b)}. (84) 

The last three terms constitute a so-called dissipative term (or dissipator) and we 
will concentrate on its form. First, we use expressions ([35]) . ([37]) for operators A a (Q) 
and 4(0). Then we consider three matrix elements. Necessary computations in 
the basis of eigenstates of free hamiltonian Ha are simple though a bit tedious, in 
some cases a suitable changes of summation indices is necessary. The results of these 
calculations are as follows 

(a\Ap(n) p(t) 4(0)1 6) = 

= y^S^mg - 0) 5(io nb -Q)(a\Aj 3 \m}(n\Al l \b} p mn (t), (85a) 

m,n 

(o|4(fi) 4(0) p(t)\b} = 

= ^2 6 ^ an ~ ^ S ( u rnn -Sl){a\A^\n){n\A p \m) p mb {t), (85b) 

m,n 



(a\p(t) 4(0) Af,(Q)\b) = 

= ^5(Ld mn - 0) 5(io bn - 0)(m | 4 I n){n\Ap \ b) p am (t). (85c) 

m,n 

The computed matrix elements are plugged into equation (|84p and summation over 
frequency O is performed. After some regrouping we find that 

^ Pab(t) = -L ^^|r a/3 (w ma ) 5(a; n6 - cj ma )(a | | m)( b \ A a \ n)* p mn (t) 

a,/3 rn,n 

- \ T al3 (uJan) S(u m n - u a n){n\ Ap \ m ) ( n \ A a \ a )* p m b(t) 

- \ T a p(uj mn ) 5(us bn - u mn ){ri\ Ap\b){n\A a \m}* p a m{t)}. (86) 

Going further, we use the evenness of Kronecker delta in the first term, while the 
presence of the deltas in the second and third terms allows us to change arguments 
in the elements of matrix T a p. Next, we denote 

K(am,bn) = ^ ^ T a p(uj ma ){ a \ Ap \ m ) { b \ A a \ n )*, (87) 
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Due to these, we rewrite formula (|86p as 

^ Pab{t) ^ = ^2 Si^ma ~ ^nb) K(am,bn) Pmn(t) 



i ^2 5{io mn - LOan) K(nm,na) p m b(t) 



2 

m,n 



~ 1 ^ W6 " K(nb,nm) p am (t). 

m,n 

Let us note the specific symmetry of this expression. Further analysis depends on 
whether the eigenfrequencies of the hamiltonian Ha are degenerate or not. We also 
note that Kronecker deltas in the second and third terms are correspondingly given 
as 8(oj mn — uJan) = S(uj ma ) and 5(uJbn — w mn ) = 5(u>bm), which allows one to perform 
summation over n. However, one has to be careful because eigenfrequencies uj n can 
be degenerate. 

5.5 Degenerate eigenfrequencies 

To account for the possible degeneracies, let us write the hamiltonian of the consid- 
ered system A in the following form 

9n 



H A = ^2huj N ^2\Nn){Nn\, 



N n=l 

where N is the main quantum number which distinguishes energy levels (energy 
multiplets), while n = 1,2, ... , gjy, are subsidiary quantum numbers. Is is obvious 
that ujn 7^ for N ^ M. Certainly, the nondegenerate case follows immediately 
and it corresponds to g^ = 1, then subsidiary quantum numbers are unnecessary 

and can be simply suppressed. 

In the degenerate case single indices appearing in equation (f88l) must be replaced 
by corresponding pairs, for example a — ► Aa. Equation (|88p is now rewritten as 

^ PAaBb{t)\ d = X] ^ UMA ~ 0JNB ) K(AaMm,BbNn) p M mNn(t) 

Mm Nn 



\ S(comn ~ uan) K(NnMm, NnAa) pMmBb(t) 

Mm Nn 

6 ( UBN - WMJV ) K(NnBb, NnMm) p Aa Mm(t). (90) 



2 

Mm Nn 



As already noted, one immediately sees that 5{<jJmn — wan) = 8{ujma) = &MA an d 
similarly 8(lobn — &mn) = &MB-, where the last deltas are the simple Kronecker 
ones. The sum over M in the second term is trivial. We put M = A and we "land 
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within multiplet A", hence we change m = a". Analogously, in the second term 
M = B and m = b" . Therefore, we have 

^ PAaBb{t) = ^2 ^2 S&MA ~ Unb) K (AaMm, BbNn) pMmNnit) 



Mm Nn 

~ 2 



1 Y,Y; K( - NnAa "' NnA ^ PAa"Bb{t) 
Nn a" 

I E K(NnBb, NnBb") p AaBh „ (t). (91) 



2 

Nn b" 

In two last terms matrix elements do not depend on quantum numbers Nn, hence 
we can denote 

n(Aa, Bb) = ^ K(NnAa, NnBb). (92) 

Nn 

This allows us to write equation (|92p in the form 

4f. PAaBb{t)\ d = ^2 E 3(u>ma ~ ojnb) K (AaMm, BbNn) p Mm Nn(t) 

Mm Nn 

- I J2 K ( Aa "> Aa ) PAa»Bb(t) ~ ±J2 K ( Bb > Bb ") PAaBV'it). (93) 
a" b" 

Let us consider this equation in some more detail. First, we take A = B (and 
correspondingly b —* a'). This yields the equation of motion for "quasi-population" 
- matrix elements taken within just one energy multiplet. Then, the first term in 
right-hand side contains 5{uma — una) = S(ujmn) = $mn- The sum over N = M 
is trivial (n — ► m!) and we have 

PAaAa'(t)\ = ^ K(AaMm, Aa! Mm') pMmMm' (*) 



Mmm' 



dt 

~ ^^2n(Aa",Aa) p Aa »Aa>(t) - \ J2n(Aa', Aa") p Aa Aa"(t). (94) 



This equation connects "quasi-populations" with other ones. The first sum contains 

the term with A = M and this term represents elastic (energy conserving) processes. 

The remaining terms (with M ^ A) corresponding to nonelastic transitions. In this 

case, the environment B serves as a reservoir which gives or absorbs the energy. The 

terms in the second line describe the "escape" from multiplet A to other ones. 
To discuss coherences, we assume A ^ B, which implies to a ^ ^b- The Kronecker 

delta in ([93]) can be rewritten as 5{ojma — ^na) = 5(uimn — ^ab)- Since ujab 7^ 0, 
we also get ujmn ^ 0. If we assume that all energy distances are different (that is 
ujab 7^ umn for different pairs A, B ^ M,N) the considered delta can give unity 
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only when A = M and B = N (which entails m — > a' and n — > b'). Then, Eq. 
reduces to 



d_ 
dt 



PA*Bb(t)\ = ^2Y1 K(AaAa',BbBb') p A a'Bb'(t) 

a' b' 

~ ±J2< Aa "' Aa ) PAa"Bb{t) ~ \ J2< Bb ' Bb ") PAaBb"(t). (95) 



So the coherences between two multiplets A ^ B couple only with coherences from 

just these multiplets. 

Obviously for the nondegenerate case "small" indices play no role - they can 

be suppressed. Then, instead of equation (|9"3j) for "quasi-populations" we get an 

equation for genuine populations 



d_ 
dt 



PAA{t) 



J2K(AM,AM) pMM(t) - K(A,A) P AA(t). 



(96) 



M 



Similarly for coherences, Eq. (|95p yields 
d 



dt 



PAs(t) 



K(AA, BB) PAB (t) 



k(A,A) + k(B,B) 



p AB (t). (97) 



These examples indicate that ME for matrix elements of the reduced density operator 
possess quite a specific symmetry which probably can be further investigated. This, 
however, goes beyond the scope of the present work. 



6 Auxiliary sections 

6.1 Preservation of normalization 

Any density operator, so also the reduced one for subsystem A must be normalized, 
that is, we require that Tr A {p A (t)} = 1. This has a simple consequence 

|Tr A W(t)}=Tr,{^}=0. (98) 

Clearly the hamiltonian part (the commutator) preserves the trace, which follows 
from cyclic property. Hence we must check the second - dissipative part of our ME. 
One may ask at which stage of our derivation such a check should be made. In 
principle, this can be done at any stage. In this section we shall do so twice. Once 
for standard form (|83p . and for ME (|93p in the energy basis. 
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Standard form 

Taking ME in its standard form (|83p we need to compute the following trace 



n a,f3 

- ±r a p(n)At(n)Ap(n) PA (t) 

- § r a p(n) PA (t) Ai(n) Ap(n)] }, (99) 

and show that it vanishes, ie., ts = 0. The trace is a linear operation, so then 



ts = EE [ r ^( n ) T ^{Ap(n) PA (t) 4(n)} 

- ± r af3 (n) Tr a {aUq) Ap(n) PA (t)} 

- \Y ap {n)Tr A {p A {t) ^(0)^(0)}]. (100) 

Cyclic property allows one to see that all three traces are equal. Therefore, tg = 
and we conclude that preservation of the normalization for ME in the standard form 
is proved. 

ME in energy basis 

In this case we check the trace preservation for Eq. (|94p . with a' = a. We need to 
compute 



ts = E E K(AaMm,AaMm') p MmMm ,(t) 

Aa Mm,m' 

- lJ2J2 K ( Aa "> ACL ) PAa"Aa{t) ~ \ k{A(1 , Aa" ) p AaAa n {t) . (101) 



Aa a" Aa 



In the first term we use definition of the parameter k (see (|92p ). In the two next ones 
we notice that indices a and a" concern the same multiplet A, so the summation 
range is also the same. We can interchange a — > a" and obtain 



ts = K ( Mm i Mm ') PMmMm'(t) 

Mm,m' 

PAaAa" 

(*) " 5 ^2^2<Aa,Aa") 

PAaAa" 

(t). (102) 



Aa a" Aa 



The second and third terms are identical and cancel out with the first one (names 
of summation indices are irrelevant). We have shown that in the energetic basis the 
trace of the reduced density operator for subsystem A is preserved. In other words, 
the derived ME preserves normalization. 
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6.2 Hermiticity of the reduced density operator 

The next necessary property of any density operator is its hermiticity. If the equa- 
tion of motion for p A (t) is identical with the similar equation for pA(t), then the 
same equations must yield the same solutions, this means that p A (t) = p A (t)- Free 
evolution is given by the hamiltonian term (—i/K) [H A + -HlSi PA{t) ] which poses 
no problems due to the commutator properties. One needs to investigate the dissi- 
pative part of ME. As in the previous section we perform such a check for ME in 
standard form and for the one in energy basis. 

Standard form 

We take the hermitian conjugate of the dissipative part of ME 

ir*^) [At(n)Aa(n), PA (t)} + }, (103) 



f! a,/3 

~ 2 



because conjugate anticommutator is equal to the anticommutator of conjugated 
operators. We know (see ()75|) ) that matrix T a p is hermitian. Interchanging the 
indices a <-> @ we get 

Jt ^ (t)) L. = + ^ E E{ P \{t) A a (n) 

il>(0) [AMsi)A p (Si), P A (t)} + }, (104) 



fl a,f3 



We see that the equation of motion for p A is formally identical with standard form 
(f83|) of ME. Thus, p A = pa - hermiticity is preserved. 



ME in energetic basis 

Hermiticity of the density operator means that PAaBb = P*BbAa- ^ entails, that the 
equation of motion for the element p%bAa mus t be the same as for PAaBb- Starting 
from ME ([93]) we look for a corresponding equation for p* BbAa - First we need to 
change the indices in (|93p (remembering that corresponding changes must be made 
for summation indices in all terms) and then we perform complex conjugation. In 
this manner we find 

J t PBbAai t )\ d = Yl Yl 6 ( U MB ~ una) K* (BbMm, AaNn) p* MmNn {t) 

Mm Nn 

- i Y,K*(Bb",Bb)p* Bb „ Aa (t) - \ Y,K*(Aa,Aa")p BbAa „(t). (105) 

b" a" 
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Next, we need to consider the conjugated quantities K* and k*. By definition (|87|) 

5(^MA ~ uJNB)K(AaMm, BbNn) = 

= ^Y, T ^ UMA ^ Aa \ A p\ Mm )( Bb \ A <*\ iVn >*- ( 106 ) 

We take complex conjugates, use hermiticity of matrix r Q( g and we notice that the 
presence of the Kronecker delta allows the change of the argument in T. Interchang- 
ing the summation indices a <-> (5 we have 

5(uj M A - cj NB )K*(AaMm, BbNn) = 

= ^J2 T ^ UNB ^ Bb \ A p\ Nn ^ Aa \ Aa \ Mm y- ( 10? ) 

Comparing this relation with definition (|87p we see that 

S(LJ MA -u NB )K*(AaMm,BbNn) = d(u) MA -u NB )K(BbNn,AaMm). (108) 

Next, we deal with the parameter K*(Aa", Aa). In the above relation we substitute 
Aa -> Nn, Mm -> Aa" , Bb -> iVn and A^n -> Aa. Then 

<5(a;A;v-w J 4Ar)-^*(iV'nAa",iVnAo) = 8{uan -UAN)K{NnAa,NnAa"). (109) 

Obviously Kronecker deltas are equal to one, so they are unimportant. Using this 
result in the definition (|92j) of the parametr k we get 

K*(Aa" , Aa) = ^ K* (NnAa" , NnAa) 

Nn 

= ^ K (NnAa, NnAa") = n(Aa, Aa"). (110) 

Nn 

Returning to the analysis of formula (|105|) . we use the proven relations (|108j) and 
(jllOl) . At the same time, in the first term in the right hand side we interchange the 
summation indices Mm «-> Nn. Moreover, we recall that Kronecker delta is even. 
Thus, we have 

loma - ^nb) K(AaMm,BbNn) p NnMm (t) 

Mm Nn 

- \ Y,<Bb,Bb")p Bb „ Aa {t) - \ Y,<Aa",Aa)p BbAa „(t). (Ill) 

b" a" 

Comparing this result with Eq. ([93|) we find that when we replace PAaBb by p* B bAa 
(consequently in all the terms) then we will arrive at (jllip . This equation is formally 
identical with Eq. (|93|) . hence p = p\ what we intended to show. Our ME preserves 
hermiticity of the reduced density operator of subsystem A. 
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6.3 Reduced density operator in interaction picture 

In this section we will show that the reduction of the density operator is invariant 
with respect to the choice of the picture. The definition ([6]) of the reduced density 
operator and ([9]) for the transformation to the interaction picture imply that the 
operator pA{t) (in the Schrodinger picture) is expressed as 

PA (t) = Tv B {p AB (t)} = TT B [e- lH ^ h Q AB {t)e iHot / h }. (112) 
We note that the free evolution of each of the subsystems written as 

I <PA(P) ) ® I ^b(O) ) — > I <p A (t) ) ® | ip B (t) ), (113) 
can be expressed with the aid of the operator 

exp (-^H t\ = exp (-~H A t \ <g> exp (-^H B t\ , (114) 

because both hamiltonians are fully independent and commute. In Eq. ([112)1 we 
compute the trace only over reservoir variables, so we can write 

PA (t) = e ~ iH At/h Tib | e -iH B t/h ~ ab ^ el H B t/h j e iH A t/h_ (n5) 
Cyclic property of the trace yields 

e iH A t/H pA{f) e - lHA t/H = ^ { ~ AB{t) } (n6) 

Left hand side represents the reduced density operator in the interaction picture (it 
depends solely on the variables of the subsystem A). Hence, we have 

g A (t) = Tv B {g AB (t)}. (117) 

This is formally identical with the definition of the reduced density operator in the 
Schrodinger picture. The relation between the reduced density and the total one 
is the same in both pictures. In other words, reduction of the operatora p A (t) = 
Tt_b { PAB^P) } is invariant with respect to the change of the pictures. 

6.4 Existence of two time scales. Discussion of approximations 

In this section we follow and (hopefully) try to to elucidate the discussion presented 
by Cohen- Tannoudji et al [6]. The previously given remarks apply also here, so we 
again refrain from frequent citations. 

Order of magnitude of time T A 

The key role in our considerations is played by the assumption ()17p . This is the 
requirement 

t b < At < T A . (118) 



S.K, J.C-K 



29 



ME - Tutorial 



6 AUXILIARY SECTIONS 



30 



In other words we assume that there exist two, quite distinct, time scales. Firstly, 
let us try to estimate the time T A which characterizes the evolution of system A due 
to the interaction with reservoir. To find such an estimate we use Eq. (|32p . that is 

/ l \2 pt+At rh 

&g A (t) =l — \J dti J dt 2 Tv B 



V AB (h), [V AB (h), gA(t)®ff B ]] J, (119) 



where we employed the introduced properties of the reservoir. We also recall that 
the main contribution to the integrals comes from a thin belt (of width tb lying 
below the diagonal t\ = t 2 , see Fig. [TJand its discussion). This allows us to estimate 
the integrand as follows 



Tr_e 



V AB (tt), [V AB (t 2 ), ~ QaT±b{v 2 o b ) = V 2 q a . (120) 



Hence, left hand side of Eq. (|119p is estimated by 

AT ~ ¥ TB v eA > (121) 

because the area of the integration region is estimated by the product T#Ai (area 
of the belt under the diagonal t\ = t 2 ). Introduced parameter V characterizes the 
"strength" of the interaction between the reservoir and system A. The factor which 
multiplies g A in (|121|) has (according to (|16|) ) the sense of the inverse of time T A . 
Therefore, we obtain an estimate 

T A ~— • or r ^~v^' (122) 

Condition for existence of two time scales 

What is the condition of the existence of two time scales? The estimate of T A given 
in (|122p allows us to find such a condition. Let us look upon condition t b <C T a 
more carefully and introduce the estimate (|122l) . This yields 

ti 2 Vt b , s 

t b « WFb - (123) 

The last inequality is the sought condition of existence of two time scales. If we 
denote £1 A b = V/h, then we can write Q A bt~b « 1. So the interaction must be 
characterized by such Bohr frequency £1 A b that during the time interval of magni- 
tude of tb its influence on system A is negligibly small. 
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Justification of weak coupling approximation 

We already mentioned (see the discussion of Eq. ([14")) ). that it is possible to iterate 
von Neumann equation - accounting for higher order corrections would increase 
accuracy. We can estimate these higher order terms in the same manner as done 
just above. For example, for the third order term we have 



At 



~ JpT 2 B Q A , (124) 



because times t\, t% and t% must be close to each other (with accuracy of the order 
of rg). Then the region of integration has volume of the order of r B At. Due to 
Eq.([l22]) we get 



A Q/ 



At 



(3) Vt b 1 _ 1 _ > \ 

~—Tw u "Ti>*< (125) 



where the last inequality follows from (|123|) . The obtained estimate shows that the 
third order iteration (and similarly higher ones) are indeed negligible. Obviously this 
holds provided the condition (|123|) is true and ensures the existence of two distinct 
time scales. 

Neglecting p core i (Born approximation) 

Moving from Eq.p4|) to (|32p we have neglected initial correlations between systems 
A and B. These correlations built up at earlier moments t' < t. This corresponds 
to the assumption that at some earlier moment to (to < t) both systems were 
uncorrelated. This happens, for example, when the interaction was switched on at 
an instant tQ. So the correlations described by Qcorei need time t — to to appear. If 
the correlations exist (g core i 7^ 0) then averaging of the term linear in interaction 
(as in expression (|3ip) would not give zero. Hence, Qcorei 7^ would result in the 
appearance of the linear term. Moreover, this would also automatically modify the 
quadratic term in (133p . Let us estimate the magnitude of this modification (which 
is due to earlier interaction) 



Aq a 1/1 



At At \ ih 



\ 2 ,-t pt+At _ 

/ dh dt 2 (VAB(t 1 )V AB (t 2 )}B. (126) 

/ J —00 J t 



The integrand contains correlation functions of the reservoir. Hence the integrand 
would be practically zero for \t\ — t%\ > tb- The integration runs effectively from 
t — tb to t + tb- Therefore, using condition (|122p . we estimate 

&QA l_ Yl 2 V^TB TB _ J_ TB , , 

At ~ At' tf TB h 2 'At T A ' At' 1 ' 

as the integrals are nonzero within the interval of the length of the order of tb- 
If tb *C At then the correction is small (main contribution to the evolution of q a 



S.K, J.C-K 



31 



ME - Tutorial 



6 AUXILIARY SECTIONS 



32 



is of the order of I/T4, which is quite larger). The key assumption that tb <C At 
allows us to conclude that the correlations between system A and B which were built 
before moment t do not significantly change the evolution of QA(t), their influence 
is restricted to the moments from a very short interval (t, t — tb)- New correlations, 
within a much longer interval (t,t + At), are building up and have an effect on the 
evolution of 5a (t). Initial correlations have small significance and hence it is justified 
to neglect them. 

Discussion of the secular approximation 

Secular approximation consists in replacing the function F(Q ' — f2) (defined in (|66p) 
in ([68]) by Kronecker delta, which leads to Eq. ([70l) . Our discussion of this replace- 
ment does not rise any doubts when |0' — 0| <C (At) -1 , because then F(Vl' — 17) is 
practically unity. On the other hand for \Q' — Q\ 3> (At) -1 the function F(Q' — O) 
is practically zero. The only problem is to justify the neglecting of the terms for 

which \Q' - Q\ ~ (At) -1 . 

To explain this point, first use the fact that the free evolution of matrix elements 

p a b is governed by 

J t PaV = ~ i<a|[W]|&> = -&w£?, (128) 
where the small Lamb shift (due to His) is ignored. The solution is simple 

p(f (t) = e — Vif(O). (129) 

Now, we analyze the dissipative term which is given by (|88|) . We do not discuss 
the nuances connected with possible degeneracies. Obviously, we can write 



SbkPmk and p am = S ^5 a kPkm, (130) 



Pmb 



which we use in (188 jl . We also interchange indices k «-> n in the second term and 
similarly, in the third one we first interchange m <-> n and then k <-» to. The result 
is 



-77 Pab(t) = y^5(uj ma -UJ nb )K(am,bn)p mn (t) 
at * ' 



m,n 



y~] $bn K^mk ~ u ak ) K(km, ka) p mn {t) 

\ ^ 5 am 5(uj bk - uj nk ) K(kb,kn) p mn (t). (131) 



1 
2 

k,m,n 



2 

k,m,n 



We note that 5 bn implies b = n, and then uj nb = 0. Therefore 

^bn S(iO m k — LO a k) = ^bn <5(<^ma) = <^bn ^(w m a ~ ^nb)i (132) 
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since uj n b = and changes nothing. Similarly we have 

5am 5(uibk — U n k) = 5 am 5(— LO n b) = 5b n ^{^>ma ~ ^nb)- (133) 

Finally, we note that 6(uj ma — w n &) = 6(10^ — u> mn ), because Kronecker delta is even. 
After these manipulations Eq. (|13ip can be written as 



d_ 

dt 



Pab(t) = ^2 ^( Wab ~~ ^mn)^K(am, bn) - \ Sbn ^ K(km, ka) 



^{kb, kn )} Pmn(t). (134) 

k 

The expression in braces is denoted as M am ^ n and we have 
d 



dt 



Pab{t) = y~] S(u ab - U mn )M am < bn p mn (t). (135) 



m,n 



This is a specific form of ME, but useful in the discussion. However, in the degenerate 

case some care must be exercised and renewed considerations might be necessary. 

We proceed to the discussion of conditions and / or possibilities of neglecting the 
terms for which |O y — f2| ~ (At) -1 . Interaction with the reservoir certainly modifies 
the free evolution of p\* (t) = e~ tuJab p ab (0). If Bohr frequencies of the oscillating 
elements p a {, i p mn are such that \u> a b — u> mn \ 3> I/T4, then mutual couplings between 
these elements are quickly averaged to zero (interfere destructively) before time Ta 
elapses and before the interaction has enough time to affect the evolution. In other 
words, if \u} a i)—u mn \ differs much from 1/Ta then the coupling between corresponding 
matrix elements will have small (weak) effect. This is the situation similar to the 
one encountered in perturbation theory. Namely, when the energies \E a — Eb\ 3> 

V a b = (a I V I b) then the perturbation has small (usually negligible) effect. 

Since, by assumption T4 3> At the discussed situation corresponds, in fact, to 

the relation \uj a b — uJ c d\ ~ (Ai) . This, in turn means, that such terms have little 

influence on the evolution of the operator pA(t)- Such terms are neglected while 

passing from Eq. (f68l) to (i70j) . Thus the last of our approximations is justified. 



6.5 Vab — Vab ~ nonhermiticity of operators A a i X a 

In our considerations we have adopted the interaction hamiltonian between the 
system A and reservoir B in the form Vab = A a <g> X a , where operators A a and 
X a do not have to be hermitian. Certainly, the full interaction hamiltonian must be 
hermitian, so we conclude that it must contain operators A a , X a and their hermitian 
conjugates Aa, X a . Constructing linear combinations we can always transform the 
interaction hamiltonian into Vab = Y2a ® -^a > where the primed operators are 
hermitian. 

We shall illustrate this with a simple example. Let the interaction hamiltonian 
be of the form 

Vab = A®X^ + At ® X, (136) 
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where operators A and X are nonhermitian, while the full hamiltonian is clearly 
hermitian. We define new operators 

Q = ^( A + A ^)> Q = ^( x + X ^)> 

P = -L(A-Ai), p = 7f( x - xt )' (137) 

which are evidently hermitian. Expressing operators A, X and their conjugates via 
the new ones, we obtain 

Vab = \{q-ip)®{Q + %P) + \{q + w)®{Q- %P) 

= q®Q + p®P. (138) 

This interaction hamiltonian is expressed as a sum of products of hermitian opera- 
tors. Hence construction of the interaction hamiltonian with nonhermitian operators 
is allowed. One can always build necessary combinations. However, in some practi- 
cal applications it is much more convenient to use nonhermitian operators than the 
linear combinations. 

6.6 Vanishing average (X a ) B 

In the main part of the lecture we assumed that Eq. (|29p holds, that is the average 
( X a )b = Tr^ { X a p B (t) } = 0. We have stated that it is not really restrictive. We 
will show that it is true. This is so, because we can always shift the energy scale. 
To see this, let us write 

Vab = Yl Aa ® ( Xa ~ ( Xa 

a 

= Y A a ®X a - (X a ) B (A a ®l B ), (139) 

a a 

where (X a ) B = TtbIctb-Xq,} is a number not necessarily equal to zero. Then we 
have 

(Va B )b = A 4( X *)b ~ (X a ) B )=0, (140) 

a 

which holds no matter whether numbers ( X a ) B are zeroes or not. Full hamiltonian 
can then be written as 

H A b = H A ® 1b + 1a® H B + V AB 

= H A ®1 B + U® H B + V AB + Y (X a ) B (A a ®l B ) 

a 

= [H A + Yl (X a ) B A a ]®l B + 1 A ® H B + Vab- (141) 
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Rescaled interaction term (the last one) has zero average (as in (|140p ). This is 
achieved by the redefinition of the energy scale in system A - via redefinition of the 
hamiltonian Ha- We conclude that the assumption that the averages (|29p vanish is 
not really restrictive, but simplifies the computation. 

6.7 Commutators of operators A a (Q) 

In the main sections we have introduced the operators A a (Sl) defined by relation 
(|35p . The hamiltonian of system A is of the form Ha = Yin n)(n\. It is not 
difficult to find the commutator [Ha, A a (Q)]. Directly from the definitions we 
obtain 

[H A , A a (U)] =[J2 

hw n \n){n\, ^S^ba - SI) | a)(a\ A a \ b)(b\\ 

n a,b 

= y^jug - u b ) 5(u ba - 0)| a)(a\A a \ b){b\ 

a.b 

= - hSl^25{uj ba -Sl)\a)(a\A a \b)(b\ = -hQA a (Q), (142) 

a.b 

which ends the calculation. Conjugation changes sign, so that 

[H A , Al(Sl)] = kSlAKSl). (143) 
Heisenberg equation of motion follows from formula (I142p . and it is 

ihj t A^\Sl) = [A£\Sl\H A ] = hSlAW(Sl). 

After integration we obtain Aa(SY) = e^A^Sl) which agrees with (|44p 
we present one more relation 

[H A , Ai(Q)A (3 (Q)]=AUsi)[H A , A p {Sl)] 

+ [H A , Al(n)]A p (Sl) = 0, (145) 

which follows immediately from the derived results. 

6.8 Explicit form of correlation functions G a p(r) 

Correlation function of the reservoir was defined in (|47p or (|53p . By assumption, 
reservoir hamiltonian Hb and the corresponding density operator o B commute, so 
they have a common set of complete and orthonormal eigenstates | z ) . Let us cal- 
culate the trace in (f53l) in the chosen basis 

G a p(r)=Tr B {xUr)X p a B } = Tr B { e «W« *t e~ iH ^/ h X p a B } 

= J2( z I eiHBT ' h X l eT lHBTth \i){i\X^a B \z). (146) 
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In Eg . ([22]) we denoted the eigenvalues of ob by p(z), hence 

Gap (r ) = £ p(z) z\Xl\£)(Z\X \z), (147) 

*,€ 

with cj 2 = E z /h , and tJ z g = w z — 

Expression (|147|) shows that the correlation function G a p(r) is a complicated 

superposition of functions which oscillate with Bohr frequencies u) z e. Reservoir is 

assumed to be large, the discussed frequencies are densely space (quasi-continuous). 

If time r is sufficiently large the oscillations interfere destructively (average out to 

zero). We can expect that reservoir correlation function decay quickly when time 

r = ti — t% increases. Characteristic decay time is denoted by tb and assumed to 

be, by far, the shortest time characterizing the system A + B. When r > tb the 

correlation may be neglected. 



7 Summary 

In this summary we describe practical steps needed in the construction of the ME 

for specified physical systems. 

The first step consists in precise definition of the system A and of the reservoir 

B. We need to specify their free hamiltonians Ha and Hb and (at least sometimes) 

their eigenenergies and eigenstates. Then we define the interaction hamiltonian in 

the form 

v AB = Y, A °® x * = J2 A i® x l ( 148 ) 

a a 

where A a , X a are (correspondingly) operators of system A and reservoir. We stress 
that these operator need not be (separately) hermitian. It suffices that the full 
interaction hamiltonian is hermitian. We also need to specify the density operator 
gb describing the state of the reservoir. It is worth remembering that operator Hb 
and ob commute. This implies that the reservoir is in the stationary state. In the 
second step of ME construction we build (identify) the following operators 

A a (n) = S("ba -n)\a)(a\A a \b){b\. (149) 

a,b 

The following matrix elements are computed in the third step 

/*oo /*oo 

W aP (Q) = / dr e iUr G aP {r) = / dr e inr Th?{Xt (r) X p a B ). (150) 
Jo Jo 

They are seen to be partial Fourier transform of the reservoir correlation functions. 
Reservoir operators are taken in the interaction picture 

X a (t) = e iHBt/h X a e - iHBt/h . (151) 
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Coefficients W a p{VL) are then employed to construct two hermitian matrices 

r afS (n) = w aP (n) + w$ a (n), 

A a(S (Q) = 1 [ W af3 (Q) - W$ a (n) ] . (152) 

We note that r a ^(f2) is a positive-definite matrix and can be computed directly as 
Fourier transform 

/oo ^ poo 

dr e iQr Tr B {xi(T)Xpa B } = / dr e iQT G a p(r). (153) 
-oo J-oo 

Parameters r a/ g(Jl), in practical applications, are more important than A a p(fl,). 
Explanation will be given later. The separate expression for elements A a p(Q) is 

-i r roo /*oo 

A a p(n) = -[J dr Tr B {xi(r)Xpa B } - J dr Tr B {xlx^r) a B } 

(154) 

Hence, calculation of coefficients W a p(VL) can be usually omitted. 

Final construction of the proper ME is the fourth and the last step. The above 
given quantities allow us to write the ME as 

j t PA(t) =-^[H A + H LS , p A (t)] 

+ p£E r a/»( n ){ WM*)4(n) 

n a,/3 

- i[4^)^(J2), PA (t)] + }, (155) 
where the so-called Lamb-shift hamiltonian Hls is given as 

H ™ = \ E E A Q/3 (f2)4(f2)^(n). (156) 

Energy shifts of the system A which are due to the presence of Hls i n the hamilto- 
nian part, are usually quite small and frequently negligible. This explains why the 
role of matrix A a/ 3 is usually less important than that of matrix r^. 



S.K, J.C-K 



37 



ME - Tutorial 



REFERENCES 



38 



References 

[1] C. Cohen- Tannoudji, B. Diu, F. Laloe, Quantum Mechanics, 
Wiley-Interscience, New York 1991. 

[2] R. Alicki and K. Lendi, Quantum Dynamical Semigroups and Applications, 
Lect. Notes Phys 717, (Springer, Berlin Heidelberg 2007). 

[3] H-P. Breuer, F. Petruccione, The theory of open quantum systems, 
Oxford University Press, 2002. 

[4] K. Hornberger, Introduction to decoherence theory, |arXiv:quant-ph/06121 18v2. 

[5] J. Preskill, Lecture notes on quantum computation, 
|http://www.theory.caltech.edu/^preskill/ph229| 

[6] C. Cohen- Tannoudji, J. Dupont-Roc, G. Grynberg, Atom-photon interactions, 
Wiley, New York 1992. 



S.K, J.C-K 



38 



